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Abstract 

Collaborative optimization is a new design 
architecture specifically created for large-scale 
distributed-analysis applications. In this approach, 
a problem is decomposed into a user-defined number 
of subspace optimization problems that are driven 
towards interdisciplinary compatibility and the ap- 
propriate solution by a system-level coordination 
process. This decentralized design strategy allows 
domain-specific issues to be accommodated by dis- 
ciplinary analysts, while requiring interdisciplinary 
decisions to be reached by consensus. The present 
investigation focuses on application of t.hfc. collabo- 
rative optimization architecture to the multidisci- 
plinary design of a single-stage-to-orbit launch ve- 
hicle. Vehicle design, trajectory, and cost, issues are 
directly modeled. Posed to suit the collaborative ar- 
chitecture, the design problem is characterized by 
95 design variables and 16 constraints. Numerous 
collaborative solutions are obtained. Comparison of 
these solutions demonstrates the influence which an 
a priori ascent-abort, criterion has on development, 
cost.. Similarly, objective-function selection is dis- 
cussed, demonstrating the difference between mini- 
mum weight, and minimum cost, concepts. The oper- 
ational advantages of the collaborative optimization 
architecture in a. multidisciplinary design environ- 
ment. are also discussed. 
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Figure 1 . Industry-team Reusable Launch Vehicle 
concepts. 

Introduction 

For several years, NASA has been examining 
various Earth-to-orbit transportation options with 
the goal of reducing operating costs relative to the 
current- tl,;S. launch fleet. [1, 2, 3, 4]. Many of these 
solutions have focused on fully reusable systems em- 
ploying various levels of advanced technology [5, 6]. 
Although a. wide range of options have been ex- 
amined, including single and t.wo-st.a.ge systems us- 
ing rocket, and/or air-breathing propulsion, current, 
emphasis has been placed on single-stage-to-orbit, 
rocket, powered vehicles [7, 8, 9, 10]. At. present., 
several single-st.a.ge-t.o-orbit. concepts are being stud- 
ied by industry-led teams with support, from NASA 
in the Reusable Launch Vehicle (RLV)/X-33 pro- 
gram [10, 11]. Some preliminary design concepts 
from this study are illustrated in Fig. 1. The design 
of a. single-st.a.ge-t.o-orbit. launch vehicle is a. com- 
plex, multidisciplinary process which is character- 
ized by thousands of design variables and nonlinear 
constraints. A complete design requires analysis of 
aerodynamics, propulsion, weights and sizing, per- 
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formance, heating, controls, operations, cost, and 
others [9]. While it is vital that each of these aspects 
be addressed at the conceptual level, the use of a 
design architecture that provides rapid performance 
of this multidisciplinary analysis while allowing the 
analyses to evolve as design maturity increases is 
equally imperative. 

Enabled by the computational advances of the 
past few decades, a suite of multidisciplinary analy- 
sis and optimization architectures have emerged [12, 
13]. Collaborative optimization is a new design ar- 
chitecture whose characteristics are well-suited to 
large-scale, distributed design [14]. The fundamen- 
tal concept behind the development of this archi- 
tecture is the belief that disciplinary experts should 
contribute to the design decision process while not 
having to fully address local changes imposed by 
other groups of the system. To facilitate this decen- 
tralized design approach, a problem is decomposed 
into subproblems along domain-specific boundaries. 
Through subspace optimization, each group is given 
control over its own set of local design variables and 
is charged with satisfying its own domain-specific 
constraints. Communication requirements are mini- 
mal since knowledge of the other groups’ constraints 
or design variables is not required. The objective 
of each subproblem is to reach agreement with the 
other groups upon values of the interdisciplinary 
variables. A system-level optimizer is employed to 
orchestrate this interdisciplinary compatibility pro- 
cess while minimizing the overall objective. This 
decomposition strategy allows for the use of exist- 
ing disciplinary analyses without major modification 
and is also well-suited to parallel execution across a 
network of heterogeneous computers [15]. 

Development of the collaborative optimization 
architecture is described in Refs. [14, 16, 17]. Com- 
parison to other optimization architectures is made 
in Refs. [14, 17]. The architecture has been success- 
fully used to solve several analytic test problems [14, 
18], trajectory optimization problems [14, 16], and 
aircraft design problems [15, 19]. The present inves- 
tigation focuses on the application of this architec- 
ture to the multidisciplinary design of a single-stage- 
to-orbit launch vehicle. Additional insight into this 
application is provided in Ref. [14]. 

Nomenclature 

A nozzle area, ft 2 
c subspace constraints 

F system-level objective 
Wing normal force, N 
g system-level constraints 


Isp 

specific impulse, s 

lh 2 

liquid hydrogen 

LHC 

liquid hydrocarbon 

LOX 

liquid oxygen 

RP 

kerosene 

Sref 

wing area, ft 2 

T 

thrust, N 

T/W 

thrust-to-weight ratio 

X 

interdisciplinary design vars 

X 

domain-specific design vars. 

y 

interdisciplinary outputs 

z 

system-level design vars. 

a 

angle of attack, deg 

AY 

ideal velocity 

subscripts 

c 

computed 

e 

exit 

si 

sea-level 

vac 

vacuum 

superscripts 

* 

subspace solution 

** 

system-level solution 


Disciplinary Analyses 

In this analysis, design of a single-stage-to-orbit 
launch vehicle includes specification of the ascent 
trajectory, as well as determination of the subsystem 
weights and sizes, and assessment of the technology 
maturation and development costs. From the results 
of a previous study, an aerodynamically viable shape 
is modeled [20]. Propulsion system characteristics 
to be identified include the liftoff thrust-to-weight 
ratio, nozzle area ratio and fuel-to-oxidizer mixture 
ratios. Two mixture ratios require specification as 
the propulsion system is operated in different modes. 
From liftoff, liquid hydrogen and liquid hydrocarbon 
are both burned as fuel with liquid oxygen; while 
during a later portion of the ascent, only the hy- 
drogen/oxygen mixture is used. The increased bulk 
density of such a dual-fuel concept has been shown to 
provide significant dry weight reductions [21]. The 
time to transition from mode 1 propulsion to mode 
2 is also optimally determined. Development cost is 
the minimization variable which, in the present anal- 
ysis, is defined to include the vehicle hardware design 
and development cost incurred under a full-scale de- 
velopment program. Other cost elements such as 
program management, fees, reserves, and software 
are not included. 
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Historically, preliminary launch vehicle design 
has not included cost, issues within the design loop. 
Instead, the pertinent engineering analyses are inte- 
grated and optimized to find either a minimum dry 
weight or minimum gross liftoff weight vehicle. A 
costing analysis is subsequently provided and, time 
permitting, some iteration may be performed. Be- 
cause the goal of the RLV program is to provide 
technology development and demonstration of a low- 
cost., reliable, space transportation system [11], this 
traditional approach is no longer appropriate. In 
this case, it. is imperative t.o consider cost, issues at. 
the preliminary design stage [22, 23]. Hence, in this 
investigation, the impact, of including cost, within the 
preliminary design process is discussed for a. repre- 
sentative, single-st.a.ge-t.o-orbit. concept.. 

The vehicle is sized to deliver and return a. 
25,000 lb payload to the Space Station following 
launch from the Eastern Test. Range at. the Kennedy 
Space Center. For this analysis, the Space Station 
is assumed to be in a. 220 n.rni. altitude orbit, with a. 
51.6 degree inclination. The single-st.a.ge-t.o-orbit. ve- 
hicle is flown into a. 50 x 100 n.rni. altitude orbit, with 
the correct, inclination; on-board propellant, is used 
to transfer to and circularize at. 220 n.rni. Burnout, 
constraints on altitude, flight.-pa.t.h angle, and incli- 
nation are enforced as are maximum inflight, nor- 
mal force, angle of attack, pitch rate, and dynamic 
pressure limits. An extension limit, is placed on the 
dual-position rocket, nozzle. In addition, technology 
maturation costs are constrained along a. represen- 
tative funding profile [24]. 

Propulsion Analysis 

Propulsion system pa.ramet.rics were supplied by 
Pratt and Whitney based on characteristics of the 
Russian RD-701 dual-fuel engine [21]. This system 
can burn either a. hydrogen/kerosene mixture (mode 
1) or pure hydrogen (mode 2) as fuel. During mode 
1, hydrogen is included in the fuel mixture to pro- 
vide nozzle cooling and increased Isp. To further 
enhance performance, the dual-fuel enginfeis fitted 
with a. dual-position nozzle. After a. regression anal- 
ysis, this parametric data, can be used as shown in 
Fig. 2. Given the two nozzle area, ratios and fuel- 
t.o-oxidizer mixture ratios, numerous engine param- 
eters are computed. These parameters include sea- 
level engine t.hrust.-t.o- weight., specific impulse, and 
maximum allowable thrust, (required inputs to the 
weights and sizing analysis) as well as vacuum thrust, 
and nozzle exit. area, (required trajectory inputs). A 
2:1 exit. area, limit, is placed on the allowable ex- 
tension of the dual-position nozzle to accommodate 
packaging of multiple engines on the vehicle base. 


r— Nozzle area ratios 
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|— Fuel/oxidizer mixture ratio, 
mode 2 
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Propulsion 
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constraint 
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Propellant bulk density 


Figure 2. Propulsion analysis inputs and outputs. 


GLOW, vehicle T/W| | iftoffi 
landed weight, base diameter 
Initial launch direction, Sref 
Pitch angles 

T„ac. V l s Pvac settings 
Propulsion/engine transition Mach numbers 
% LH2, % RP, model 
Mixture ratio, mode 2 


— Other weights & sizing inputs 

• Aero coefficients 

• Atmospheric model 

• Mission information (orbit, launch site, flight 
duration, normal-force boundary duration) 
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Figure 3. Trajectory analysis inputs and outputs. 


Trajectory Analysis 

To analyze the ascent, flight.-pa.t.h, a. three 
degree-of-freedom trajectory analysis is performed 
with the Program to Optimize Simulated Trajec- 
tories (POST) [25]. Within POST, the equations 
of motion are numerically integrated from an ini- 
tial to a. terminal set. of state conditions. Within 
the present, investigation, the vehicle is treated as a. 
point-mass, Earth rotation and obla.t.eness are mod- 
eled, and the 1976 standard atmosphere (no winds) 
is used. As shown in Fig. 3, the required set of 
trajectory inputs includes vehicle (e.g., gross liftoff 
weight., liftoff t.hrust.-t.o-weight. ratio, and aerody- 
namic coefficients and reference area.) as well as tra- 
jectory parameters (pitch angle history, launch az- 
imuth, and the propulsion-syst.em transition Mach 
number). This domain-specific analysis is respon- 
sible for evaluation of the inflight, and terminal 
constraints, computation of the vehicle mass ratio, 
and determination of the required fuel and oxidizer 
masses (weights and sizing inputs). Terminal con- 
straints on altitude, velocity, and flight.-pa.t.h angle 
as well as maximum inflight, dynamic pressure, an- 
gle of attack, pitch rate, and normal force (wing siz- 
ing constraint, based on landed weight.) limits are 
enforced. 
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Fuel/oxidizer mass fractions 
% LH2, % RP, model 
Desired mass ratio 
Desired landed wing loading 

i— Other weights & sizing inputs 

• Vehicle geometry 

• Propellant densities 

• Mission information 
(crew size, duration, orbit) 


^ Weights & Sizing J 

r 

Subsystem 

weights 

Mass ratio GLOW c 

wing loading Sref c 

constraints Base diameter c 

Figure 4. Weights and sizing analysis inputs and 
outputs. 

r- Technology readiness levels 
Subsystem weights 



Technology Development 

maturation cost 

cost 

constraint 

Figure 5. Cost analysis inputs and outputs. 

Weights and Sizing Analysis 

The; Configuration Sizing program (CONSIZ) 
developed at NASA Langley Research Center is used 
to size the vehicle and determine subsystem weights. 
This sizing process is performed to meet vehicle 
mass ratio and landed wing loading constraints. As 
shown in Fig. 4, the liftoff thrust-to- weight, fuel- 
t.o-oxidizer mixture ratio and fuel/oxidizer mass, as 
well as several propulsion system parameters are 
required input to CONSIZ. With the exception of 
the liftoff thrust-to-weight and the fuel-t.o-oxidizer 
mixture ratio, all of these inputs are computed by 
one of the other two disciplinary analyses. In ad- 
dition to dry weight, CONSIZ computes the gross 
liftoff weight, reference aerodynamic surface area 
and landed weight (each of which is a required tra- 
jectory input). Numerous subsystem weights are 
also required as input to the cost, analysis. 


Cost Analysis 

As depicted in Fig. 5, the cost, model used in 
the present, investigation consists of two sets of cosi- 
est, ima.t.ing relationships. The first, set., technology 
maturation costs, provide an estimate of the invest- 
ment. required to advance a. given subsystem from 
it.0 present, technology readiness level to a. higher 
technology readiness level within five years (the pre- 
sumed beginning of the development, cycle). As pre- 
sented in Fig. 6, the NASA technology readiness 
scale was used as a. measure of subsystem maturity. 
Data, from Ref. [24] was used to compile these re- 
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Level 5 
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Level 8 
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Technology concept and/or application formulated 
Analytical and experimntal critical function and/or 
characteristic proof-of concept 
Component and/or breadboard validation in 
laboratory environment 
Component and/or breadboard validation in 
relevant environment 

System/subsystem model or prototype demonstration 
in a relevant environment (ground or space) 

System prototype demonstration in a space 
environment 

Actual system completed and "flight qualified” 
through test and demonstration (ground or space) 
Actual system "flight proven" through successful 
mission operations 


Figure 6. NASA technology readiness levels. 


la.t.ions. For this estimation, the vehicle subsystems 
are partitioned into nine segments: avionics, struc- 
tures, electromechanical actuation, electrical conver- 
sion and distribution, auxiliary propulsion, prime 
power, main propulsion, propellant, tanks, and ther- 
mal protection. The total technology maturation 
cost, is defined as the sum of the cost, of maturing 
these nine elements. 

The second set. of nonlinear, vehicle-specific 
cost.-est.ima.t.ing relationships were developed con- 
sistent. with the subsystem breakdown structure of 
CONSIZ. These parametric relations predict, subsys- 
tem development, cost, as function of weight., com- 
plexity, and technology readiness level. Additional 
relationships for subsystems which do not. rely on 
advanced technology (e.g., landing gear) were in- 
cluded to obtain an estimate of the total vehicle 
development, cost, assuming a. full-scale development, 
program. The development, cost, relations were for- 
mulated through application of the Parametric Re- 
view of Information for Costing and Evaluation- 
Hardware (PRICE-H) multivariate model [26]. This 
model has been used to generate integrated cost, and 
schedule estimates for numerous industry and gov- 
ernment. aerospace programs. The model provides a. 
computational method for deriving cost, estimates of 
electronic and mechanical hardware assemblies and 
systems. 


Interdisciplinary Coupling 

From Fig. 2 to Fig. 5, it. is clear that solu- 
tion of this problem requires an iterative approach 
as the each analysis each requires inputs which are 
computed by another discipline. For example, one 
must, ensure that, the reference aerodynamic sur- 
face area, resulting from the vehicle sizing process 
(Sref e of Fig. 4) is the same as the reference aerody- 
namic surface area, used to compute the aerodynamic 
forces (Sref of Fig. 3). Interdisciplinary compati- 
bility must, also be achieved in regard to the gross 
liftoff weight., landed weight., base diameter, mass 
ratio, fuel/oxidizer mass, and each of the propulsion 
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Interdisciplinary 

Variable 

Prop. 

Traj. 

Weights 

Cost 

Gross liftoff weight 


input 

output 

Aero ref. area 


input 

output 

Landed weight 


input 

output 

Base diameter 


input 

output 

Vehicle liftoff T/W 


input 

output 

% LH 2 , mode 1 

input 

input 

input 

% RP, mode 1 

output 

input 

input 

Mixture ratio, mode 2 

input 

input 


Vac. thrust, mode la 

output 

input 


Vac. Isp, mode la 

output 

input 


Vac. thrust, mode lb 

output 

input 


Vac. Isp, mode lb 

output 

input 


Vac. thrust, mode 2 

output 

input 


Vac. Isp, mode 2 

output 

input 


Noz. area, retracted (a) 

output 

input 


Noz. area, extended (b) 

output 

input 


Mass ratio 


output 

input 

LH 2 mass fraction 


output 

input 

RP mass fraction 


output 

input 

Engine sea-level T / W 

output 


input 

Sea-level Isp, mode 1 

output 


input 

Allowable thrust ratio 

output 


input 


Table 1 . Interdisciplinary coupling present in 
launch vehicle design problem. 

discipline outputs. The liftoff thrust-to-weight and 
fuel-to-oxidizer mixture ratios must be treated in a 
similar fashion since these parameters are input to 
more than one analysis. The interdisciplinary cou- 
pling structure of this launch vehicle design problem 
is shown in Table 1. Here, the input /output struc- 
ture of the propulsion (Prop.), trajectory (Traj.), 
and weights, sizing, and cost subspaces is depicted. 
In this table, the weights, sizing, and cost analy- 
ses are grouped together since they are integrated 
within the same subspace in the collaborative archi- 
tecture solutions. 

Potential Optimization Architectures 

Numerous optimization architectures have been 
proposed for launch vehicle design. In this section, 
some of the potential solution strategies are dis- 
cussed. The design problem, as posed within the 
collaborative optimization architecture, is then pre- 
sented. 

The disciplinary tools used in this analysis were 
originally developed as stand-alone programs, each 
operated by a disciplinary expert. To obtain a cred- 
ible vehicle, a design team was required to manually 
iterate among these disciplinary analyses. In many 
cases, “optimization” was performed through trade- 


studies in which the parameters were varied one at 
a time [27, 28]. 

Improvement over this one-variable-at-a-time 
approach has been achieved using response surface 
methods (RSM) [9, 20, 21, 29, 30], In a RSM 
strategy, feasible designs are computed at numer- 
ous points in the design space and a surface is fit to 
these points. Optimization is then performed on this 
approximate representation of the design space. Al- 
though use of RSM is a significant improvement over 
one-variable-at-a-time trade studies, the method suf- 
fers from several drawbacks including the require- 
ment to produce numerous feasible design candi- 
dates and the limitation to problems characterized 
by a relatively small number of variables. For larger 
applications, as of interest in this work, the use of 
RSM is not a viable option. 

Numerous gradient-based optimization archi- 
tectures are possible for solution of this launch ve- 
hicle design problem. For example, use of system 
sensitivity analyses is presented in Ref. [31] for a 
similar problem (cost issues not addressed). Within 
the launch vehicle design community, the traditional 
multidisciplinary design optimization approach is to 
integrate the appropriate disciplinary analyses in a 
nonhierarchic iterative loop with a single, system- 
level optimizer [32, 33, 34, 35, 36]. In this standard 
optimization approach, interdisciplinary compatibil- 
ity is enforced through some form of loop conver- 
gence criteria. As a result, numerous analysis evalu- 
ations are required to produce a single design candi- 
date. An all-at-once architecture has also been pro- 
posed for launch vehicle design [23, 36]. In this ap- 
proach, the iterative loop of the standard approach 
is removed through the use of auxiliary variables and 
compatibility constraints [15, 37]. As demonstrated 
in Ref. [36], use of a sequential-analysis, all-at-once 
optimization architecture may significantly reduce 
the required number of function evaluations. In a 
subsequent section, the standard and all-at-once so- 
lution strategies presented in Ref. [36] are compared 
to the use of collaborative optimization. 

Collaborative Optimization 

Figure 7 shows the launch vehicle design prob- 
lem as posed to suit the collaborative optimization 
architecture. Here, the design is decomposed into 
three subspace coordinated by a system-level opti- 
mization procedure. Note that the weights, sizing, 
and cost issues are all accommodated within a single 
subspace. This decomposition strategy was selected 
because the weight and cost models are tightly cou- 
pled through numerous subsystem weights. For ex- 
ample, in the present model, development cost is a 
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Figure 7. Collaborative optimization architecture 
for launch vehicle design. 
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Figure 8. Interdisciplinary negotiation within the 
collaborative optimization architecture. 


Optimization 

Level 

Design 

Var. 

Nonlinear 

Const.. 

System 

23 

3 

Trajectory subspa.ee 

47 

9 

Propulsion subspa.ee 

4 

1 

Weights, sizing, cost. sub. 

21 

3 


Table 2. Collaborative optimization launch vehicle 
design characteristics. 

function of 17 subsystem weights. Hence, in this 
case, analysis integration is useful in reducing the 
degree of interdisciplinary coupling. 

Posed to suit the collaborative archifjgcture, the 
problem is characterized by 95 design variables and 
16 constraints. These design variables and con- 
straints are partitioned among the system-level and 
subspaces as listed in Table 2. Further insight into 
the subspace and system-level problems is provided 
in Ref. [14]. 

Numerous collaborative solutions were obtained 
for this launch vehicle design problem. In each case, 
the subspace bounds for the interdisciplinary vari- 
ables were not selected consistently. For example, in 
the trajectory subspace, the vehicle liftoff T/W was 
required to be between 1.0 and 1.3; whereas, a range 
of 1.0- 1.5 was placed on this interdisciplinary input 
within the weights, sizing, and cost, subspace. Simi- 
larly, the percentage of LHo in mode 1 had bounds 
of 0-15, 3-20, and 3-15 within the various subspaces. 
These bound inconsistencies were included to sim- 
ulate the architecture’s application in a team envi- 
ronment where the various disciplinary specialists 
are each responsible for a particular domain-specific 
subproblem. 

Prior to optimization, the subspace and system- 
level problems were scaled such that the design vari- 
ables, constraints, and objective function were of 


order one. However, in this investigation, no at- 
tempt was made to provide a twice-continuously dif- 
ferentiable model. A majority of the sources which 
contribute nonsmoothness are within the trajectory 
subspaee-t.he ascent pitch profile which is modeled 
by discrete control points connected by linear seg- 
ments as well as tabular linear interpolation of the 
aerodynamic and atmospheric properties (1976 stan- 
dard model). 

As an example of the collaborative convergence 
process, Fig. 8 depicts the interdisciplinary negotia- 
tion that occurs between the trajectory subspace, 
weights-sizing-cost. subspace, and system-level on 
the appropriate value of the wing aerodynamic ref- 
erence area (Sref). This figure shows the domain- 
specific convergence histories for 3 system-level iter- 
ations. In the first, iteration, the system-level pro- 
poses an Sref value of 4000 ft 2 . After meeting all of 
its domain-specific constraints, the trajectory sub- 
space returns with a. request, to change this value 
to 3950 ft. 2 . Similarly, the weight.s-sizing-cost. sub- 
space returns with hopes of increasing this value to 
4100 ft 2 . Based on gradient, information provided 
at. the solution of each subspa.ee, a. system-level step 
is taken which alters the system-level value Sref to 
3955 ft. 2 (and the value of the other 22 system- 
level variables). In this case, the subspa.ces are 
able to remain disciplinary constraint, feasible while 
providing better agreement, on this value of Sref. 
In the third iteration, the value of Sref is reduced 
even further with consent, from the two domain- 
specific analyses. Through repeated collaboration, 
the system-level optimizer orchestrates the interdis- 
ciplinary compatibility process. Note that in Fig. 8, 
cold-start, subspa.ces were used for illustrative pur- 
poses. Hence, within each system-level iteration, 
subspa.ee optimization was initialized from a. fixed, 
domain-specific point, (an Sref of 4000 ft. 2 in the 
trajectory subspa.ee and an Sref of 5000 ft. 2 in the 
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Figure 9. Collaborative optimization system-level 
convergence profile. Minimum development cost, 
concept. 


weights-sizing-cost subspace). In the results which 
follow, subspace optimization is initialized from the 
prior solution (warm-start.) to reduce the computa- 
tional requirements of repeated subspace optimiza- 
tion. 


Application of the Collaborative 
Optimization Architecture 

In this section, optimal solutions of several ver- 
sions of the basic launch vehicle design problem are 
presented. Each case was obtained through appli- 
cation of the collaborative optimization architec- 
ture. To obtain these solutions, the refined archi- 
tecture of Ref. [14] was used. This implementation 
includes use of the linear system-level objective re- 
finement. as well as the warm-start and slack- variable 
subspace refinements presented in Ref. [14]. The 
system-level Jacobian is obtained from the use of 
post-optimality information estimated at subspace 
solutions as discussed in Refs. [14, 38]. The sequen- 
tial quadratic programming algorithm, NPSOL [39], 
was used t.o provide system-level and subspace opti- 
mization. CPU times quoted are based on use of a. 
Silicon Graphics Challenge L machine outfitted with 
R8000, 90 MHz processors. 

Minimum Development Cost Concept 

The collaborative optimization convergence his- 
tory for this solution is presented in Fig. 9. The 89 
system-level iterations shown required 181 subspace 
optimizations of each of the three domain-specific 
analyses. The solution obtained compares favorably 
(within 0.1% in the system-level objective function) 
with the solution presented in Ref. [23] where an 
a.ll-a.t.-once optimization strategy was used. As dis- 
cussed in Ref. [14], this convergence profile, where 
the system-level objective initially drops and then 
approaches the solution from below as system-level 
feasibility is achieved, is similar to collaborative so- 
lutions obtained for other problems. 
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Figure 10. Collaborative optimization convergence 
history: system-level values of vehicle mass ratio and 
wing aerodynamic reference area.. 



Figure 11. Collaborative optimization convergence 
history: system-level values of main engine sea-level 
T/W and vehicle gross liftoff weight. 

For the solution shown in Fig. 9, 181 sets of 
subspace optimization werp performed. Initialized 
from a. fixed initial guess ( cold-start. ), a. call to the 
subspaces requires approximately 1-3 hours. At this 
rate, a. cold-start, solution requires approximately 1-3 
weeks. Using warm-start subspaces (restarting from 
the previous domain-specific solution with knowl- 
edge of the prior optimum active set, Lagrange mul- 
tipliers, and Hessian of the La.gra.ngia.n), the solution 
presented in Fig. 9 required approximately 4.5 days 
of computer time. Hence, warm-starting the sub- 
spaces provides a. dramatic advantage. This level of 
efficiency gain is possible because the subspaces are 
tasked with solving a. related sequence of subprob- 
lems. 

Figures 10, 11, and 12 demonstrate the conver- 
gence behavior of six of the interdisciplinary vari- 
ables. These figures depict, the system-level values of 
the vehicle mass ratio, wing aerodynamic reference 
area., main engine sea-level T/W, vehicle gross liftoff 
weight, vehicle liftoff T/W, and the mode 2 mixture 
ratio. As shown, the system-level variables tend to 
exhibit either the same behavior as the system-level 
objective or a. more damped oscillation about, the 
optimal value. 

The flight path portion of this optimal solution 
is illustrated in Figs. 13 and 14. The ascent, begins 


7 





Figure 15. Minimum development-cost, vehicle de- 
scription, no low-altitude ascent, abort constraint. 


Figure 12. Collaborative optimization conver- 
gence history: system-level values of vehicle liftoff 
T/W and oxidizer/fuel mixture ratio during mode 2 
propulsion. 



Time, s 

Figure 13. Optimal flight profile (altitude, velocity, 
and acceleration). 



Time, s 

Figure 14. Optimal pitch-angle and inflight con- 
straint. profiles. 


with a. 400 ft. vertical rise to clear the launch facil- 
ity. Because the vehicle is characterized by a. low 
liftoff T/W (1.036), this rise takes approximately 20 
s. This vertical flight, segment, is followed by a. max- 
imum pitch-rate segment, in which the vehicle tries 
to attain its maximum-lift, orientation. The pitch- 
rate is limited t.o 5 deg/s t.o reflect, control issues 
which are not. modeled in this analysis. As shown in 
Fig. 14, during this segment, of flight., angle of attack 
(a) increases until it. reaches the allowed maximum 
of 20 deg. Flying at. this a, the normal force builds 
until a. limiting load is reached. The vehicle rides 
this normal-force boundary through peak dynamic 
pressure which for the optimum flight, path is about. 
990 psf. Hence, at. this solution, the dynamic pres- 
sure limit, of 1000 psf is not. active. During this phase 
of flight, at. approximately 10 kft. , the back pressure 
losses are low enough that, the dual-position nozzle 
is extended t.o gain propulsive efficiency. As the noz- 
zle extension is performed, the vehicle acceleration 
initially decreases as a. result, of flight, through the 
transonic regime. At. about. 50 kft., as the dynamic- 
pressure decreases, the vehicle comes off the normal- 
force boundary but. continues t.o accelerate towards 
3 g’s. The vehicle reaches 3 g’s while in the dual- 
fuel propulsive mode and uses engine throttling to 
maintain this level of acceleration. In this analysis, 
transition of all seven engines from a. dual-fuel to 
single-fuel mode is performed instantaneously. This 
instantaneous change in thrust, results in the large 
decrease in acceleration shown in Fig. 13 at, roughly 
225 sec and an altitude of about, 175 kft,. Operating 
in a, single-fuel mode (LHo), the vehicle accelerates 
back to 3 g’s and holds this acceleration level un- 
til reaching orbit,. Note that, if the nozzle transition 
been performed sequentially, a, slightly better result, 
could have been achieved. 

This minimum development, cost, configuration 
is sketched in Fig. 15. As shown in this figure, the 
vehicle is characterized by a, dry weight, of 1.831 x 
10 5 lb, a, gross liftoff weight, of 2.139 x 10 b lb and a, 
length of 181.6 ft,. Note that, the optimum vehicle 
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Figure 16. Collaborative optimization system-level 
convergence profile. Minimum development cost, 
concept (vehicle liftoff T/W > 1.2). 

T/W is 1.036 at liftoff. This T/W level, which is 
significantly lower than the 1.2 value typically ob- 
served, is t.h#;subject. of the following section. 

How Many Engines? 

With the system-level vehicle liftoff T/W lower 
bound increased to 1.2, the collaborative optimiza- 
tion architecture was re-run, giving the system- 
level performance profiles shown in Fig. 16. Here, 
86 system-level iterations are required to obtain 
system-level convergence. Comparison of the solu- 
tions presented in Figs. 9 and 16 demonstrates that 
requiring a vehicle liftoff T/W >1.2 induces a 4.5% 
increase in development cost.. 

Generally, a. lower bound on the vehicle liftoff 
T/W is included at. the preliminary design stage to 
provide sufficient, engine-out. capability during the 
first, moments of ascent, (vehicle T/W increases with 
time). In this case, 1.2 was chosen because the ve- 
hicle was designed with 7 engines. Hence, with loss 
of a. single engine at. liftoff, the vehicle would still 
have propulsive control a.ut.horit.y-t.he T/W would 
become 1.03 [40]. Inclusion of this engine-out. capa- 
bility is the primary reason that, many launch vehi- 
cle design concepts have a. liftoff T/W of 1.2 (e.g., 
Refs. [7, 9, 20, 21, 31, 36]). Unfortunately, this in- 
creased flexibility does not. come without, a. price, as 
the minimum development, cost, concept, would prefer 
a. lower liftoff T/W. 

To further complicate this issue, the number of 
engines is a. significant, driver on vehicle cost, and re- 
liability. As the number of engines increases, flight, 
reliability concerns also increase. To maintain a. cer- 
tain level of system reliability, e.g., 99%,. bach en- 
gine in a. 7 engine cluster must, have a. reliability of 
99.86%. With just. 3 engines, this individual engine 
reliability is reduced to 99.66%. Since reliability is a. 
significant, cost, driver [1, 6], low-cost, concepts such 
as those envisioned within the RLV program may not. 
have the luxury of including this engine-out. abort. 


constraint.. For example, with 3 engines, incorpora- 
tion of a. liftoff engine-out. capability induces a. se- 
vere penalty since, the lower bound on the allowable 
liftoff T/W would be increased to 1.5. In this case, 
this liftoff engine-out. requirement, may have to be 
discarded. Although this is the philosophy selected 
in the single-st.a.ge-t.o-orbit. designs of the present, in- 
vestigation, it. is clear that, this issue requires further 
study. 

Objective Selection in Launch Vehicle Design 

The concept, of a. “best.” or optimal config- 
uration has different, meanings in different, design 
groups. In launch vehicle design, minimum weight, 
concepts have traditionally been sought.. In some 
cases, liftoff gross weight, is selected as the minimiza- 
tion variable [27, 28]. However, with the realization 
that, propellant, is relatively inexpensive, a. major- 
ity of the n cen i concepts have been design to min- 
imum dry weight, (e.g., Refs. [7, 9, 20, 21, 31, 36]). 
Often, these minimum dry weight, studies include a. 
claim such as, ! f Since vehicle development costs tend 
to vary as a function of dry weight, this minimum 
dry weight vehicle may be considered a minimum de- 
velopment cost concept”. However, as demonstrated 
in this section, such an assertion is not. rigorously 
true even when a. weight-based cost, model is used. 

Another objective function which has been sug- 
gested within the launch vehicle design community 
is minimum AV. This objective is sometimes cho- 
sen based on application of the “rocket, equation” . 
However, as shown in this section, minimum AV 
concepts are vastly different, from either minimum 
development, cost, or minimum weight, designs. In 
this case, it. is the use of analysis approximation (i.e., 
use of the rocket, equation instead of a. higher-fidelity 
trajectory simulation program) which leads to such 
a. dramatically different, concept.. 

Distinctions among these four different, objec- 
tives (gross weight., dry weight., development, cost., 
and AV) are examined through application of the 
collaborative architecture. In each of these designs, 
the vehicle liftoff T/W is allowed to vary in the range 
1.0-1. 5. The minimum development, cost, solution, 
described in the previous section, is used to normal- 
ize the other optimal results. 

Design characteristics of these four optimal con- 
cepts are listed in Table 3. As shown by the fourth 
column of this table, at. the vehicle level, the mini- 
mum AV concept, stands apart, from the other 3 de- 
signs. This system is over 50%) more expensive than 
the minimum development, cost, concept, and more 
than 25%) heavier than the minimum dry weight, case 
(while yielding only a. 4-6%) improvement, in AV). 
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Normalized 

Design 

Characteristics 

Dev. 

Cost. 

Objective Function 
Dry Liftoff 

Weight. Weight. 

AV 

Dev. cost. 

1.0 

1.001 

1.040 

1.510 

Dry weight. 

1.0 

0.987 

1.035 

1.262 

Liftoff weight. 

1.0 

1.007 

0.925 

1.196 

AV 

1.0 

0.999 

0.987 

0.944 

Body length 

1.0 

0.993 

1.030 

1.052 

Wing area. 

1.0 

0.993 

1.043 

1.239 

Liftoff T/W 

1.0 

1.004 

1.082 

1.448 

Mass Ratio 

1.0 

1.008 

0.893 

0.959 


Table 3. Optimal launch vehicle comparison. 



development dry gross liftoff A V 

cost weight weight 

Figure 17. Optimal launch vehicle comparison- 
main propulsion system. 

Furthermore, the vehicle’s wing area is more than 
20%) larger than each of the other wings. With the 
exception of this concept, the other three concepts 
share a top-level similarity, the strongest correlation 
being between the minimum development cost, and 
minimum dry weight configurations. However, even 
among these first, three concepts, the maximum vari- 
ance in dry weight, is 3.5%), the maximum variance 
in development, cost, is 4.0%), and the maximum vari- 
ance in gross liftoff weight, is 6.7%). 

The optimal liftoff T/W and mass ratio for each 
of the four concepts is also listed in Table 3. From 
the liftoff T/W comparison, it. is clear why the mini- 
mum AV concept, is so different.. In this case, the op- 
timization process attempts to minimize the ascent, 
losses (e.g., gravitational, aerodynamic, nozzle back- 
pressure, and thrust- vectoring) by achieving orbit, 
as quickly as possible. As shown in Fig. 17, to ac- 
complish this, the propulsion system size is greatly 
increased and a. larger RP propellant, mass fraction 
(LHC) is required. In this case, the vehicle liftoff 
T/W is at. its upper bound of 1.5. Had this variably 
been allowed to increase even further, the vehicle 
would have continued to grow, reaching enormous 
proportions. 


Figure 17 shows that, subsystem differences are 
not. limited to the minimum AV concept.. For ex- 
ample, while a. dual-fuel system is allowed, the min- 
imum gross weight, concept, prefers LHo as its sole 
fuel. In this case, the dense RP propellant, is elim- 
inated such that. Isp is increased while the initial 
launch weight, is reduced. This disregard for bulk 
density is in stark contrast, to the minimum dry 
weight, case which relies on bulk density to provide 
dry weight, savings [21]. The goal of minimum gross 
liftoff weight, also induces a. 8.2%) increase in vehi- 
cle liftoff T/W and a. 10.6%) decrease in vehicle mass 
ratio (seq Table 3) relative to the minimum develop- 
ment. cost, system. 

Figure 17 also shows that, the minimum devel- 
opment. cost, and minimum dry weight, concepts are 
not. equivalent.. This finding is of even greater signif- 
icance when one considers that, the present, analysis 
relies on a. weight-based cost, model. For example, 
relative to the minimum dry weight, design, the min- 
imum development, cost, concept, uses 16%) less RP 
(LHC). In this case, the minimum cost, concept, seeks 
t.o reduce propulsion system mass since this subsys- 
tem is characterized by both a. high development, 
cost, per pound and technology maturation cost.. In 
essence, the minimum development, cost, case maybe 
thought, of as a. “weighted” minimum weight, case, 
where all pounds are not. created equal. 

The results shown in this section demonstrate 
the importance of appropriate objective function se- 
lection. In particular, while minimum weight, and 
minimum cost, concepts may be similar, they are 
certainly not. equivalent, (even when a. weight-based 
cost, model is used). Differences at. the subsystem 
level were shown to reflect, the goal of the optimiza- 
tion process. Therefore, if one is truly interested in 
developing minimum cost, concepts (the stated goal 
of the RLV program), cost, considerations must, be 
factored into the preliminary design process prior to 
optimization [22, 23]. 

Operational Aspects of Architecture Selection 

In Ref. [14], the computational performance 
of the collaborative optimization architecture was 
shown to improve as the size of the domain-specific 
problems are increased without, a. significant, increase 
in interdisciplinary coupling. Hence, using compu- 
tational performance as a. yardstick, application to 
large-scale, loosely-coupled problems was suggested. 
In a. multidisciplinary design environment., computa- 
tional performance is only one aspect, of architecture 
selection. Often, in such a. setting, analysis inte- 
gration and communication requirements are much 
harder to efficiently resolve. In addition, system 


10 




Optimization 

Architecture 

Function 

Evals 

Mod. 

Time 

Comm. 

Reqs. 

Std. approach, [36] 

10482 

4 mo. 

66 

Collaborative opt. 

3125-24840 

1 mo. 

23 

All-at-once, [36] 

3182 

3 mo. 

65 


Table 4. Comparison of three multidisciplinary op- 
timization strategies for launch vehicle design. 

flexibility, analysis modularity, and resource man- 
agement are also concerns. For example, if it re- 
quires a year to set up the appropriate multidisci- 
plinary analysis system, run-time is certainly incon- 
sequential. Furthermore, if a system is integrated in 
a manner such that it is difficult to modify or ex- 
tend, it is not likely to have a long life within the 
design organization. As the number and fidelity of 
the domain-specific analyses increases, these opera- 
tional concerns becomes more significant. 

Based on the results of this paper and Ref. [36] 
(where cost issues were not addressed), Table 4 
presents the author’s experience solving a similar 
launch vehicle design problem with three optimiza- 
tion architectures. This table lists the function eval- 
uations required to reach the solution from a com- 
mon initial guess, the estimated analysis modifica- 
tion time (Mod. Time), and the interdisciplinary 
communication requirements (Comm. Reqs.) for 
each of the three optimization architectures. As a 
result of the variation in subspace problem sizes (see 
Table 2), the propulsion subspace required far fewer 
function evaluations than the trajectory subspace. 
This is the cause for the range of numbers listed in 
the function evaluation column of Table 4 for the col- 
laborative architecture. This table shows that while, 
in its present form, the collaborative architecture 
is not computationally competitive with the all-at- 
once optimization architecture, it is generally com- 
petitive with the standard optimization approach. 
A similar conclusion was reached for the large-scale 
trajectory optimization problems of Ref. [14]. Be- 
cause Table 4 lists function evaluations, the potential 
parallelization speedup inherent to the collaborative 
and the all-at-once approaches is not reflected. 

Use of the collaborative architecture offers nu- 
merous other advantages. The largest of these is the 
estimated analysis modification time. Recall that 
this analysis is based on the modification of a set 
of previously stand-alone disciplinary programs. In 
such a case, the standard and the all-at-once ap- 
proaches require a much higher degree of analysis 
modification. For the standard optimization ap- 
proach, this time is spent providing integration of 


the disciplinary analyses. Since this integration was 
performed in an “engineering” sense, the resulting 
integrated system is not extremely flexible. 

In the case of the all-at-once system, signifi- 
cant time was spent preparing the analyses for opti- 
mization. In a practical design environment, many 
disciplinary analyses are simply not set-up for effi- 
cient optimization (“design-oriented”). For exam- 
ple, in its original form, the trajectory analysis was 
not well-suited for function evaluation. Instead, tra- 
jectory evaluation was explicitly coupled with opti- 
mization. While the original analysis program was 
simply inserted into the collaborative architecture, 
significant time was spent providing a trajectory 
function evaluator for use with the other approaches. 

Communications requirements of the collabo- 
rative architecture are also minimal (see Table 4), 
resulting solely from the interdisciplinary coupling 
inherent in the problem. In contrast, both the 
standard and all-at-once approaches required signif- 
icantly more communication. In the standard ap- 
proach, 36 variables were passed between the opti- 
mizer and the set of multidisciplinary analyses. The 
remaining communication requirements involved co- 
ordination among the analyses. Similarly, in the all- 
at-once strategy, 40 variables were shared between 
the optimizer and the analyses, while 25 variables 
were passed among the analyses. In contrast, the 
collaborative architecture only requires communica- 
tion of 23 variables between the system-level and 
the subspaces. This decreased level of communica- 
tion is a direct result of empowering the subspaces 
with domain-specific decision responsibility. 

In a multidisciplinary design environment, use 
of the collaborative architecture provides additional 
operational advantages as a result of its flexibility 
and modularity. Within this architecture, both anal- 
ysis and optimization may be performed distribu- 
tively. In fact, the distributed aspects of the archi- 
tecture have been demonstrated on a set of hetero- 
geneous computing platforms in Ref. [15]. Further- 
more, in collaborative optimization, the required de- 
composition into a set of coordinated subspaces ne- 
cessitates a coarse-grained modularity. This modu- 
larity allows the system to be easily extended and 
modified. Over the lifetime of a design project, this 
flexibility may be used to adjust the fidelity of the 
disciplinary models. In addition, changes to one 
domain-specific analysis may be altered without im- 
pacting the rest of the system. In contrast, the stan- 
dard optimization approach provides very little flex- 
ibility and was found to be difficult to modify, while 
the all-at-once approach only provides a distributed 
analysis capability. 
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Summary 

In this investigation, the collaborative optimiza- 
tion architecture was used to perform multidisci- 
plinary design of a dual-fuel, single-stage-to-orbit 
launch vehicle. Vehicle design, trajectory, and cost 
aspects were directly addressed. Posed to suit the 
collaborative architecture, this design problem was 
characterized by 95 design variables (23 interdisci- 
plinary) and 16 constraints. A minimum develop- 
ment cost concept was obtained which compared fa- 
vorably (within 0.1% in development cost) to a re- 
sult produced by another optimization architecture. 
The influence of an a priori ascent-abort criterion on 
development cost and proper objective-function se- 
lection was discussed. Differences were highlighted 
between the minimum cost, weight, and AV con- 
cepts. 

The operational aspects of the collaborative ar- 
chitecture in a multidisciplinary design environment 
were presented and compared with two other op- 
timization architectures. Relative to these other 
optimization strategies, the advantages of the col- 
laborative architecture include no analysis integra- 
tion requirements, the ability to use domain-specific 
analyses which already provide optimization without 
modification, inherent system flexibility and modu- 
larity, a distributed analysis and optimization capa- 
bility, and a significant reduction in communication 
requirements. These practical advantages make the 
architecture well-suited for use in a large-scale, mul- 
tidisciplinary design environment. 
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